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Abstract 


Mathematical models based on the finite element 
method of structural analysis as embodied in the 
NASTRAN computer code are widely used by the 
helicopter industry to calculate static internal loads 
and vibration of airframe structure. The internal 
loads are routinely used for sizing structural 
members. The vibration predictions are not yet 
relied on during design. NASA's Langley Re- 
search Center sponsored a program to conduct an 
application of the finite element method with 
emphasis on predicting structural vibration. The 
Army/Boeing CH-47D helicopter v.^s used as the 
modeling subject. The objective was to engender 
the needed trust in vibration predictions using 
these models and establish a body of modeling 
guides which would enable confident future pre- 
diction of airframe vibration as part of the regular 
design process. 


Introduction 


A better capability to calculate vibration of 
helicopters is a recognized industry goal. More 
reliable and accurate analysis methods and 

computer aids can lead to reduced developmental 
risk, improved ride comfort and fatigue life and 
even increased airspeeds. An important element in 
the overall vibration calculation is the finite 
element airframe model. Under a NASA Langley 

Contract, Boeing Vertol Company performed the 
program enunciated by the title of this paper, that 
is, the planning, creating and documenting of a 
NASTRAN finite element vibration model of a 
modern helicopter. Further, test requirements 
were established and a ground shake test 
performed to validate the model. An unusual 
requirement of the contract was that each major 
step of the program be presented to and critiqued 
by the Industry, 

The contract consisted of two phases with multiple 
tasks in each phase: 

• Phase i. Planning, Creating and Documenting 
A Helicopter NASTRAN Model 

Task 1 Planning 

Task 2 Modeling 

Task 3 Test Requirements 

Task 4 Industry Critique 


Presented at the Second Decennial Specialist's 
Meeting on Rotorcraft Dynamics, AHS/NASA Ames, 
Moffett Field, CA, Nov. 8, 1984. 


e Phase II: Vibration Test to Verify a 

Helicopter NASTRAN Model 

Task 1 Aircraft Ground Shake Test 
Task 2 Industry Critique 
Task 3 Report 


Functions of the Finite Element Models (FEM's) 

The forming of FEM's have become almost routine 
for new helicopter airframes. But to step back a 
moment, why are they being formed? . . . what are 
the current uses after they are formed? . . . and 
what are the future uses as the technology 
improves and the degree of correlation advances? 
"Today's" functions of the finite element model 
(FEM) static models are shown in Figure 1. They 
are commonly used to calculate fuselapa internal 
loads. What formerly was an extensive job invol- 
ving months of effort by many Stress engineers 
has been reduced to routine running of cases once 
the FEM is prepared. Then the same model “an be 
the basis for a vibration model. 
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FUNCTION TECHNICAL OECWON IMPACT 


• USE STATIC FEM MODEL TO CALCULATE 
INTERNAL LOADS 

- THE CRITICAL LOADS ON EACH AIRFRAME 
ELEMENT TOHRMIT SIZING AND STRESS 
ANALYSIS 

i MAJOR. FEM INTERNAL DESIGN LDAOS 
USED 

- REOUCES STRESS MAN LOADING FOR 
INTERNAL LOAD CALCULATIONS 

• STATIC MODEL USED AS IASIS FOR THE 
VII RATi ON MODEL 

• IASIS FOR MODEL LEADING TO VIBRATION 
CONFIGURATION DECISIONS 
- N0> COMMONLY DONE YET 


• wu STATIC M00EL TO CALCULATE 

deflections 

• OCCASIONAL USE 


• CALCULATE AIRFRAME FATIGUE LOADS 

• NONE /ET 

- FUTURE CAFAIIUTY 



Figure 1 . Function* oi Static. Finite E Ivntnt 
M odeJti 

"Future" functions of the FEM include calculation 
of airframe fatigue loads. Field problems with air- 
frames often involve cracking of skin panels or 
stiffeners from vibratory loads. Early prediction 
and correction of such problems would be a useful 
improvement to the aircraft. 

For vibration models, the categories of functions 
can be discerned for the engineering development 
of helicopters. These are (1) guiding structural 
design so as to avoid resonance with rotor exciting 
frequencies, (2) predicting flight forced vibration 
levels, and (3) supporting design of vibration con- 
trol devices. 
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There are two "Today" functions in Figure 2 which 
are in routine use. The first row is the function 
to predict and control resonances in the basic de- 
sign. Three forcing frequencies are addressed; 
1/rev, b/rev and 2b/rev. By far the most severe 
and limiting vibration occurs at b/rev. However, 
even relatively low vibration levels at 1/rev can be 
annoying. 2b/rev levels are next in importance, 
and can be significant when seeking very low 
overall vibration. FEM's are employed at the detail 
design stages to check proximity of the lower 
natural modes to 1/rev. If analysis indicates a 
proximity which is judged to be a concern, the 
procedure would be to utilize the analysis to ex- 
plore corrective structural changes and to imple- 
ment these changes in the design before 
construction. The changes could affect both the 
structural arrangement and the structural gages. 

FUNCTION TECHNICAL DECISION HWfACT 

• PREOICT NATURAL FREQUENCY PLACEMENT ■ MINOR FOR t> RE V AND ZbHEV 

AN0 credibility FOR THESE HIGHER 

MOOES POOH SO GENERALLY LITTLE 

• MOOiEV Of SIGN BEFORE ORAWINC RELEASE IMPACT 

TO ASSURE REQuiREO NATURAL FREQUENCY 
PLACEMENT 

• MODERATE FOR l/REV ANO AE ROE LASTIC 
MOOES 

SOME RECOGNITION IN AVOIDING THESE 
PROBLEMS BEFORE DRAWING RELEASE 


• IOENTIFY STRUCTURAL MODEL • MODE RATf 

MODIFICATIONS^ ANTICIPATION OF SOME uSf f QR p R£ p* RINC HARDWARE 

NEED FOR IMPROVED TUNING AFTER REAQY FQR SHAKE TEST QR FLIGHT 

SHAKE TEST OR FLIGHT TEST U$T f r£ que NCY TUNING 


• PREDICT FORCED VIBRATION UNDER 
INOIVIOUAL UNIT LOADS <1 E INDIVIDUAL 
ROTOR LOAD DIRECTION CROWN PRESSURE) 

ANO • NONE 

• MODIFY DESIGN BEFORE DRAWING FUTURE FUNCTION 

RELEA>{ TO ACHIEVE MINIMUM FORC'D G S 


* PREOiCT UNTREATED FORCED VIBRATION 

ANO • NONE 

• MODIFY DESIGN BEFORE DRAWING - FUTURE FUNCTION 

RELEASE TO ACHIEVE MINIMUM fORCEC C $ 


• DETERMINE SIZE ANO EFFECTIVENESS OF . NONE 

VIBRATION TREATMENT DEVICES . 0CCAW0MAt AIRFRAME ABSORBER 

SIZING 


• PREOICT TREATED FLIGHT Vl|R*TIO* 

ANO * NONE 

• MODIFY TREATED OESIG.Y BEFORE ORLWIFG - 'UTURE FUNCTION 
RELEASE TO MEET FLIGHT MIRATION $PfCa 

Figure 1 , Function* Vibration Finite. Element 
Model* 

The FEM is almost certainly applied during detail 
design to check for proximity of any of the higher 
airframe modes with b/rev. But, based on today's 
perspective it is not so predictable what actions 
engineers would undertake preceding prototype 
fabrication if a proximity of concern should be 
indicated by the analysis. Two reasons for the 
uncertainty are present: (1) higher mode behav- 

ior of the airframe has been regarded as difficult 
to predict, and (2) due to weaknesses of the 
currently available tools to predict vibration levels 
of the coupled rotor/airframe system it could be 
difficult to reach a consensus on whether any 
predicted coincidence of a natural frequency with 
b/rev reflects a real problem. There has, 


however, been some use of FEM for estimation of 
the effectiveness of stiffening hardware in raising 
natural frequencies. This is perceived to be more 
dependable because only the delta frequency 
rather than the absolute frequency is used. 

For the future, it is expected that forced vibration 
from individual rotor vibratory loads and from 
combined rotor loads will be predicted on a routine 
basis. Not only will they be predicted, but the 
airframe design will be iterated before drawing re- 
lease to minimize forced vibration levels. 

Modeling Plan 

As a counterpoint to most modeling efforts, this 
program emphasized the planning of the modeling 
as the prime portion of the effort. All of us have 
modeled by spreading out the drawings and getting 
down to work, typically without a very clear idea 
of where we were headed. In contrast to this, the 
NASA Technical Monitor : nsisted on a well thought 
out plan of attack, accompanied by detailed pre- 
planned instructions, labeled "guides". These 
guides defined the modeling approach for each 
type of structure-frames, stringers, rotor shafts, 
etc. Even the documentation of the modeling had 
to be preplanned. A very extensive modeling plan 
report, Ref.(1) was published. The plan was re- 
viewed by other Industry representatives prior to 
undertaking the actual modeling. Another unique 
feature was that at the end of the modeling, dev- 
iations from the planned guides due to cause were 
reported . 

The objectives of the modeling plan were as 
follows: 

• Define guides for modeling, coding, docu- 
menting and demonstrating (1) stress (static) 
modeling, (2) mass modeling, and (3) vibra- 
tion modeling (by modification of the stress 
model). 

Establish the organization, schedule and re- 
sources for performing detailed finite element 
modeling of a CH-47D helicopter. 

• Identify and discuss the functions of finite 
element vibration models in the design 
process. 

• Provide for plan critique by the industry. 
Modeling Guides 

Guides for static, mass and vibration modeling 
were developed. These included 

• Node and element numbering 

• Frame, stringer, skin treatment. 

• Rotor shaft and transmission modeling. 

• Concentrated and distributed masses. 

• Changes from the static model to form 
a vibration model. 
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The aircraft was first divided into major areas for 
convenience in scheduling and tracking FEM activ- 
ities. For the CH-47D, the breakdown was as 
shown in Figure 3. 


FRAME CAPS 



3 U.'MTER fuselage 
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5 AFT FUSELAGE .. 
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Figu/ie, 3. 8*eafedowtn into M ajoK Aavu ioK Statxc 
Modeling 
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A logical grid and element numbering scheme was 
selected to permit traceback of the elements. 
Blocks of numbers were assigned to major sections 
as indicated in Figure 4. 


FiguAe. 5. Static Modeling Guide* - Fxame* 
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Figure. 4. Mode. and Element Numbexing Scheme 


Detail guides for modeling were described. Sev- 
eral typical guides are illustrated in Figures 5 and 

6a 


F-cguAe 6. Static ModtLLng Guide* - 8utkhe.acU, 
Veck*, and B 'xtt-Line. Beam 


The mass modeling procedure is summarized in 
Figure 7. Mass data for the aircraft were first 
compiled on a standard weights tape per 
Ml L-STD-4S1 or M I L-STD-1 374 . Masses were then 
divided into concentrated items and distributed 
items. Concentrated Items such as transmissions 
and engines were allocated to Individual NASTRAN 
nodes of the static model. Distributed items, 
structure, wiring etc., were allocated to frame 
stations by Boeing program W-17, and then manu- 
ally distributed to nodes at that station. 





ITEMS TO BE INCLUDED 



Figure 7, M <U6 Modeling Guide* 


The planning effort highlighted the fact that a 
good static model may serve as the vibration model 
with relatively small changes as shown in Figure 8. 



Flgu/it 5. Vibration Modeling Guide* - Change 
faom Static to Vib flatten TB I 

The vibration model used CQDMEM2 elements to in- 
clude the axial stiffening effectiveness of skin 
panels and webs, which were neglected in the 
static model by the use of CSHEAR elements. The 
logic was that under limit loads, the skins buckle 
and do not contribute much to axial stiffness. 
In the vibration case under 1g static loads, the 
skins are unbuckled and effective. 

Documentation 



LEVEL I OVERVMW 

1 PRIMARY FUSELAGE STRUCTURE 

2 STRUCTURE MOT MODE* '0 

3 NASTRAN MODEL 


LEVEL N AIRFRAME MAJOR EECTIOR 

1 AIRFRAME MAJOR SECTIONS 

2 NODE AND ELEMENT NUMBERING 
AND IDENTIFICATION SCHEME 


LEVEL III AIRFRAME BUB~BECTK>M BREAKDOWN 
I. COCKPIT 

2 FORWARD FUSELAGE 
1 CENTER FUSELAGE 
4 FUEL PODS 
8. AFT FUSELAGE 
B. AFT PYLON 


LEVEL IV AIRFRAME MOD* LINO DETAILS 
1 ROTOR SHAFT 
2. TRANSMISSION COVER 
1 FRAME 

4. BULKHEADS, OECKS. AND BUTT LINE BEAM 
1 SKIN AND STRINGER 
8. FLOOR 
7 FUEL POO 
I LANDING GEAR 
». ENGINES 


Ftgusiz 9. Fofunal Documentation Plan (M Static 
Modeling 

being modeled, its relationship to the major area 
and its corresponding NASTRAN model, and the 
rationale for modeling assumptions along with the 
details of section property computations. 

Indurtry Critique of Modeling Plan 

In an approach which is becoming more common in 
government supported research, other industry 
members participated in the program. 

Boeing, the prime contractor, was required to 
subcontract to other major helicopter manufac- 
turers, a series of review tasks. Bell, Hughes 

and Sikorsky were the participants. Upon comple- 
tion of the modeling plan, Boeing briefed the 
subcontractors at their own sites, and reviewed 
verbal and written commentary on what the others 
thought of the plan from their own background of 
experience. 

Examples of the comments were: 

- the use of substructuring via superelements 
was suggested for cost and time saving. 

- a more detailed mass model was recommended 

- stringer lumping to save complexity and cost 
was questioned. 


An important aspect was the documentation plan, 
Figure 9. Quite often, modeling and documentation 
are done on an "as I get to it" basis. In this 
program, all of the steps were preplanned. The 
documentation was planned at four levels: over- 

view, major sections, subsection breakdowns and 
modeling details. The documentation was to pro- 
vide a clear Illustration of each major area being 
modeled, a clear Illustration of particular details 


- the forward transmission cover model was too 
simplified 

- procedures for checking the model should 
have been defined, such as SPC checks, rigid 
body checks etc. 

This review procedure was repeated later for the 
test plan, and for the analytical correlation. 
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Actual Modeling Experience 


The static model was prepared by a senior stress 
engineer and a technician working from the draw- 
ings of the CH-47D. Figure 10 shows the final 
NA5TRAN model of the aircraft with the statistics 
indicated. 


RASTRAN HOD CL 
1,113 STRUCTURAL NODES 
S.7SI STRUCTURAL HIRIRTS 
MO. Of 



Vigwiz 10. CH-47V WA5TRAW StAuctwul Uedtl 



ViguAe, 11. Static Modeling o £ FemmAcf Rotox 
Shait and T turn* min ion CovtA 

The rotor shaft was represented by two CBAR ele- 
ments with node points at the bearing locations. 
A cruciform structure comprised of CBAR elements 
was used to model the transmission cove, . The 
cover model provided bearing node points to sup- 
port the rotor shaft, and node points at the air- 
frame attachments. Bending stiffness of the 
transmission cover legs was represented by the 
four legs of the cruciform model. 

Modeling details of a typical center fuselage frame 
are shown in Figure 12. 


The planned numbering system, previously pre- 
sented in Figure 4, was straightforward, easily 
applied, and required a maximum of only four 
digits for grid points and five digits for elements. 
In the case of the grid points, sequential num- 
bering was possible which facilitated checking for 
missing points in the listing. Capability was 
provided for independent modeling (except at 
interfaces) of the major airframe sections which is 
a necessary feature for rapid development of a 
model . 

There were disadvantages turned up. The loca- 
tions of nodes and alements were not obvious from 
the numbers. Only general location was implied by 
the block number. Any later revisions or addi- 
tions tended to disrupt the numbering sequence 
and patterns. A principal difficulty was the es- 
timation of the number block sizes. If sufficient 
space was not allocated, the numbering sequence 
was inter; upted. Estimating an adequate number 
of grid points was relatively simple, but estimating 
sufficient space for the elements was difficult. 
This situation could have been partially alleviated 
by coding the element types which then would 
have made the full block of numbers available for 
each element type. The system of using station 
numbers in the code Is probably the best, 
although it increases the *Ue of the identification 
numbers. 

Details of a typical subassembly static modeling 
task are illustrated by the model of the forward 
rotor shaft and transmission in Figure 11. 



figuAt 12. Static Modeling o{ Sta. 200 FAame 

The caps carried axial load only and were repre- 
sented by CON ROD'S. Average cap area was used 
between nodes where the cap was tapered. Cap 
areas were reduced for fastener holes, and local 
cap notches were Ignored. No portions of adja- 
cent skin or effective areas of webs were lumped 
with the caps. Webs carried only shear and were 
modeled with CSHEAR's. Web holes and stiffeners 
were ignored. 
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The model for the forward pylon upper buttline 
beams is shown in Figure 13. 






Tiguit 13. Static. Modeling oh Fomind Pylon 
UppeA Butt- Lint Beam 


The transmission support fitting at the top of the 
beam was designed to act as a truss and is mod- 
eled with axial CONROD's. Otherwise the model 
was like a frame in that caps were represented by 
CONROD's and webs by CSHEAR's. Stiffeners 
used only for web stability were not all mooeled 
(some were to break up panel sizes). 

Longerons, stringers and side skins were modeled 
as in Figure 14. Longerons were modeled as 
CONROD's using their actual areas. Stringers, 
because there a-e 36 of them on the cross-section, 
were lumped into 13 effective stringers (or lumped 
with longerons) to limit the size of the model. 
Skin panels were represented by CSHEAR's. 



FtguAe 14. Static Modeling oh Stungvu, 
LongtAOiU, and Stdi SklnA 


Modeling of effective skin near longerons and 
stringers as an addition to their area was one of a 
number of instances where the guides were vio- 
lated for cause. The logic was originally that the 
static (stress) model would recognize buckled skins 
occurring under design maneuver loads. Then 
with this buckled skin model, internal load dis- 
tributions would be obtained for detailed stressing 
of the elements. Locally effective areas of skin 
were to have been added to stringer areas for po- 
tentially improved accuracy. 

For the vibration model, the airframe was to have 
been treated as in 1g level flight without maneuver 
induced buckling. The original guide was written 
to remove the locally effective skin area from the 
stringers for the vibration model where the skins 
were to be fully effective. 

It was realized when the actual modeling was 
underway, that the labor of adding and then 
removing tnese small delta areas was not worth- 
while. The static model internal load distributions 
would not really be affected by these small delta 
areas. This change was the most significant of 
the deviations made from the planned guides. 

A demonstration run was made with the static 
model to determine whether the model generated 
reasonable (error free) results. Internal loads 
were calculated for a 3 g pull-up at a gross weight 
of 50000 pounds. Element forces, grid point dis- 
placements, and grid point force balances were 
examined. The static deflection plot for selected 
grid points illustrated in Figure 15 indicates 
apparently rational results. 
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Next, the model had to undergo certain modifica- 
tions from a static to a vibration model. One of 
these changes was the drag strut of the engine 
mount. The drag strut, Figure 16, is slotted and 
only acts under extreme maneuver and crash 
loads. It was included in the static model, but 
was removed from the vibration model. The 
inactive strut has a vibration purpose; it prevents 
the drag strut from adding a yaw stiffness 
increment which would have placed the engine yaw 
natural frequency on 3/rev. Further, since the 
forward yoke support fitting is significant in 
forming the stiffness of the engine mounting, this 
yoke was remodeled to provide better detail. Cap 
areas of the forging were modeled with CBAR’s 
and the webs with CQUAD2 shell elements. 



TiQu/ie 16* VlbKcution ModtUng S&uic.taAal Chang U 


The most important change to form the vibration 
model was the change of airframe skin from 
CSHEAR's to CQDMEM2 elements. The latter are 
membranes which provide both the skin shear cap- 
ability and are effective in adding bending area. 
The change is associated with the buckled versus 
unbuckled skin configurations of the static and 
vibration models discussed previously. 

Concentrated weights of the engines, trans- 
missions, and APU were initially distributed to the 
attachment points in the static model while pre- 
serving the mass and inertia of the overall air- 
craft. For the vibration model/ center of gravity 
grid points were introduced at the engines and 
transmissions and appropriate inertias used. 

A demonstration run was performed with the vibra- 
tion model. It was done in the free-free condition 
to represent an inflight situation. Emphasis was 
placed on the basic airframe structure by modeling 
an empty aircraft without fuel. This avoided the 
need for dealing with the nonlinear cargo and fuel 
isolation systems. The demonstration run 
included the calculation of natural frequencies and 
modes and forced response. Results of the natural 
frequency calculation are summarized in Table 1. 
Based on previous CH-47 modeling and test 
experience, these results were judged to be 
reasonable. The modeling process was reported in 
Ref. 3. 


Table. 1. HbAjcutLon Vm$n&ttoLtion Cau, 
UaXuAal 


MODE >0. 

FACOUINCV (H2) 

01 SCt 1 FT ION 

1 

1.31 

1ST LATERAL - ATT NYLON LATERAL 

2 

7.24 

ENGINE LATERAL YAH - OUT 07 PHASE 

3 

7.92 

1ST VERTICAL - ATT 7YL0R LONGITUDINAL 

4 

1.24 

ENGINE LATERAL YAH - IN PHASE 

1 

11.19 

2ND VERTICAL - PYLON LONGITUDINAL IN PHASE 

1 

12.19 

2N0 LATERAL . FU0 PYLON LATERAL 

7 

11. II 

3RD LATERAL - P.YLOi LATERAL IN PHASE 

1 

14.01 

ATT LANDING GEAR LATERAL - OUT OF PHASE 

9 

11.22 

UNDE7INE0 VERTICAL 

10 

17.41 

UNDE7INC0 LATERAL 

11 

19.20 

UN0EMNED LATERAL 

12 

20.71 

UNOEFtNEO VERTICAL 

11 

21 .41 

UNDEFINED VERTICAL 

II 

22.92 

UNDEFINED COUPLED VERT 1 CAL -LATERAL 

11 

24.79 

UNDEFINED COUPLED VEIT 1 CAL-LATIRAL 

0 

1 

b/NV 24/ RCV 

10 1 IS 20 i 21 N2 


Time and Cost 

A key question has long been, can an FEM be 
assembled and used in time to influence the design 
of a new helicopter airframe? 

This was estimated in great detail, as illustrated in 
Figure 17, and it appears that an initial vibration 
result can be obtained in 6 months from Contract 
Award. This is certainly timely, because primary 
structure releases are not completed until the 15th 
month. 
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The cost of the modeling is 4430 manhours or 5% of 
a typical 85000 manhour airframe design effort. Of 
this 5%, 4% is for the static model - an activity 
that is becoming routine by Stress, and only 1% 
more to obtain the first vibration model results. 
Beyond this point vibration iterations of the 
design will add to the cost, but will certainly be 
cost effective if it provides a well tuned fuselage 
prior to manufacture. 
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Test Plan 


In addition to flight vibration measurements, two 
categories of ground tests can be identified as 
a means of evaluating a finite element model of a 
helicopter airframe for vibration analysis; namely, 
static deflection tests and shake tests. The ground 
test approaches have two significant advantages: 
(1) the rotor is removed which is a great simpli- 
fication, and (2) ail applied forces can be mea- 
sured and controlled. 



Fxgu^e It. SumaJiy oi Volution Ttit Load 

Conditions 

While the deflection test was deemed to be very 
desirable, it was not performed because of cost 
limitations. 


Shake Test 


The objective of the shake test was to verify the 
NASTRAN finite element vibration model. The 
approach was to obtain detailed frequency response 
and mode shapes under conditions which exercised 
all elements of the model. These included 
excitation at both forward and aft hubs using all 
flight hub forces and moments except torque and 
covering the frequency range from 5 to 35 Hz 
(9/rev is 33 Hz). 

The planned method of excitation was to suspend 
electrodynamic shakers and the aircraft from a 
shake test gantry. Figure 19. 


VI ItT I CAL SHAHID 


INPLANE SHAH* 


Static deflection tests seem attractive because: (1) 
Inertia effects are eliminated allowing independent 
evaluation of stiffness. (2) To some extent, 
selected parts of the airframe can be loaded facil- 
it 4 ting identification of model deficiencies. On the 
nsgaiiv^ side, industry experience with complete 
airframe deflection tests is extremely limited. 
Finally, it is noted that correlation with n shake 
test directly addresses the proposed application. 

Deflecti o n Test 

The objective of the deflection test was to verify 
the stiffness modeling performed analytically. The 
approach was to obtain detailed deflection data 
under loading conditions which exercised all major 
structural elements of the airframe. These 
included bending, torsion and frame racking of the 
constant section, pylon bending, and pylon to 
constant section load path. 

The proposed deflection test loadings of Figure 18 
were designed accordingly . 



/iguAt 19. Shake Te&t Method oi Excitation 


Dual vertical shakers operating in a master/slave 
mode are driven ;n or out of ph:se to provide ei- 
ther vertical, pitch or roll excitation. In the 

vertical direction, the soft suspension of both the 

aircraft and shaker isolates the shaker from the 

aircraft except through the drive link. In the 

horizontal plane, Isolation of the shaker is pro- 
vided by the low frequency pendulum modes of the 
aircraft and shaker on the suspension cables. 
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Airframe accelerometer locations are shown in 
Figure 20. 




F^guAe 20. Shake Tut kui^tome MeatuAvnint 

Location* TjLquap. 21. CH-47P Tut Specimen in Shake Tut 

Support Fixture 


Response measurements were to be obtained at 5i 
locations in three axes. Locations correspond to 
node points of the NASTRAN vibration model. 

Pretest NASTRAN forced analysis results were to 
be compared with the shake test results. The pri- 
mary criteria for correlation was intended to be 
the forced response plots. Secondary criteria 
would be the mode shapes at the natural 

frequencies. 

Industry Critique of Test Plan 

As with the modeiing plan, an industry review of 
the Ref. 2 test plan took place. With regard to 
the desirability of the deflection test, one 
considered the cost to outweigh the benefit. Two 
pointed out that modal parameters including 
damping should not be neglected. Two noted that 
the selection of hub mass effect is an important 
aspect of the test. And two reminded us that 
rotor shaft and drive system free play may have 
a significant impact on results. 


Ground Shake Test and Correlation 


The test specimen was the second prototype of the 
VCH-47D helicopter, Figure 21, 



As per the test plan, the aircraft was suspended 
at the rotor heads in a large structural steel fix- 
ture which also supported the rotor head shakers. 
A low frequency suspension, all less than 2. Hz, 
was employed for both the iircraft and shakers. 
Three linear vibratory forces and two moment: 
were applied at each rotor head. Selection of 
force levels was based on practical considerations 
including sufficient magnitude of response, shaker 
stroke limits and stable behavior of the susoended 
shakers. Results were obtained in the form of 
transfer function plots and mooe shapes for each 
excitation, Figure 22. 


Figu/te 22. Format oi Typical Shake Tut Ruutts 


For e< excitation an extensive matrix of forced 
response plots was obtained. Figure 23 is an 
example. 

A summary of the test natural frequencies devel- 
oped from the matrix of response peaks s pre- 
sented in the Figure 24 bar chart. The shaker 
excitation which provided the best excitation is 
noted , 
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F*gu*e 23. Frequency Rtipomi SummaAy fan 
Tonwcuid Longitudinal Excitation 


From location to location, considerable scatter was 
sometimes evident in the frequency at which a 
given mode appeared. Th : c made it difficult to 
precisely define the natural frequencies. Observed 
nonlinear behavior with force level is believed to 
be at least partially responsible for the scatter in 
the peak frequencies. In the bar chart of Figure 
24, the frequency with the largest response was 
favored . 

Fc xed mode shapes for the two vertical and two 
lateral modes closest to 3/rev are displayed in 
Figures 25 through 28. All of the shapes repre- 
sent the total forced response normalized by the 
maximum deflection. The first response sh*pe at 
T1.7 Hz (Figure 25) is dominated by the longi- 
tudinal pitch motion of the forward pylon with a 
smaller in-phase motion of the aft hub. Motions of 
the two hubs are balanced by an essentially rigid 
body motion of the remainder of the aircraft. 
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FujuAe 14 . StumaAy o^ Ttit Wo tiJuxJL Etitqutnciib 


figtinz 25. FoAced Mode Shape at 11,7 Hz with. 

FoAwatd Hub Longitudinal Excitation 



EiguSLt 16, FoAced Mode Sh ape at 11,6 Hz with 
FoAxoaAd Hub LateAal Excitation 


The characteristic of the 12.6 Hz mode, Figure 26, 
is essentially that of a classical second torsion 
mode. A relatively large lateral/roll motion of the 
forward pylon is accompanied by a small in-phase 
motion of the aft pylon. The pylon motions are 
opposed by a differential lateral motion of the 
upper and lower cabin structure. A large lateral 
motion of the aft landing gear also contributes to 
the inertial balance. 
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Like the previous mode at 12.6 Hz, the response 
at 15.2 Hz, Figure 27, :? also basically a second 
torsion mode. In this case, however, the in-phase 
hub motions are opposed by wnat more nearly re- 
sembles a twisting motion of the cabin, as indi- 
cated by differential motion from le r t to right as 
well as top to bottom. Note also that the phase of 
the aft landing gear is reversed in the mode. 


shaker effects. Consequently, the total model was 
fully representative of the test configuration in- 
cluding the support fixture, the shakers and the 
aircraft and shaker suspension systems in addition 
to the basic airframe model. A differential stiff- 
ness correction was also developed and applied to 
the stiffness matr* < to include gravitational effects 
(pendulum mode^j on the suspended aircraft. 




rlguAe 28. Forced Mode Shape cut 16.2 Hz with 
fomaAd Hub VeAtlcal Excitation 


At 16.2 Hz, Figure 28, the response shape (Dis- 
played is the fundamental vertical bending mode of 
the cabin section. Bending motion of the cabin is 
opposed by large out-of-phase pitch motions of the 
pylons. 

NASTRAN Analysis of Test Configuration 

The basic airframe vibration FEM initially demon- 
strated in the free-free condition was modified to 
the test configuration. Changes to the basic air- 
frame model included incorporation of the test hub 
fixtures (hub weight and shaker beam as c embly) 
and adjustments to the mass distribution to account 
for equipment not installed. 

The total NASTRAN model incorporated several 
unique features. A persistent issue with regard 
to analytical correlation of test and analysis has 
jeen the question of the suspension system and 


Shaker and Support System Modeling 

Modeling of a typical shaker configuration is illus- 
trated in the schematic of Figure 29. The shaker 
stator mass and a portion of the ^radle assembly 
mass are located at the shaker pivot point (grid 
7011). The remaining cradle assembly weight is 
located at the cradle suspension point (grid 7012). 
The armature flexures (armature spring) connect 
the stator and the coincident armature mass. 
Motion of the armature mass is constrained to act 
along the axis of the drive rod. The drive rod, 
represented by a CONROD, is assumed to carry 
only axial loads due to the flexures oriented at 
90°. 



Figure. 29. Typical Shake K and Siupen&lon 
Modeling 
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Fiime 30. Support Uxtufiz NASTRAN Model 

The NASTRAN model of the shake test support 
fixture which weighs approximately 37,32< pounds 
is shown in Figure 30. Grid points corresponding 
to the aircraft and shaker support points are iden- 
tified. Typical modeling of the hub and shaker 
suspension which is the interface between the 
support fixture and the basic airframe is illus- 
trated in Figure 31. 
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F^gufte 3/. FotoitoAd Hub SuApeiuion Modeling 
bon VeMicaL/Pitch Excitation 


RESPONSE: FWD. HUB VERTICAL 

EXCITATION: FWD. HUB VERTICAL 



FiguAe 32. Typical Analytical R e*pon*t boK 
F*ee and Su&ptndtd Condition* 


With regard to the question of the suspension 
system and shaker effects, the support fixture is 
always likely to have modes in the test range. 
The question, therefore, can only be resolved by 
a comparison of analytical aircraft responses for 
the free and suspended conditions. Typical 
results illustrated in Figure 32 show only minor 
effects with the most significant changes in the 30 


to 35 Hz range. While these results are applicable 
only to the test equipment used in this program, 
they generally support the accepted suspension 
concept. Physically, frequency shifts and ampli- 
tude variations may result from any of the fol- 
lowing or combination of the following: 

• Coupling with shaker system 

• Minor coupling with the support fixture 

• Prestiffening of the airframe due to gravity 
preload. 

• Other coupling mechanisms in the airframe 
due to gravity preload. 

Also, it should be remarked that the theoretical 
appropriateness of representing pendulum modes 
by a differential stiffness correction, while 
plausible, has not been thoroughly explored. 

Correlation of Test and Analysis 

Conventional correlation of test and analysis for 
airframe vibration is a comparison of natural fre- 
quencies and modes first, and forced vibration 
second. In this program the criteria order was 
reversed; more emphasis was placed on the ability 
of the analysis to predict reasonable forced ampli- 
tudes throughout the airframe. Natural modes 
were in second place, although it is recognized 
that specific forced peaks and valleys follow 
natural frequency placement. If able to predict 
reasonable forced amplitudes from individual 
rotor forces, then the analysis would be a reason- 
able tool for predicting vibration arising from 
actual mixed forces and directions. 

To keep the correlation process within reasonable 
bounds, forced vibration results were presented at 
only four representative and widely separated loca- 
tions, Figure 33, each in the vertical, lateral and 
longitudinal directions. The forces for illustration 
were the forward rotor vertical, pitch and lateral 
excitations. \ single structural damping value of 
2.5% critical was assumed. 


• FOMKQ WTO* nr.TIUL, PITCH MO LATCMi [ROTATION 

• FOftcco « at four hi oar scmmteo locations 



• SUSKNKO UST CORHQURATIOH 

• MODAL DAMIRS EQUAL 2 . SI CRITICAL 


F^uAe 33. AiKbKajnt Location* and Condition* b ok 
V oiced R e*pon*e CowielatLon 
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Forced response comparisons with forward vertical 
excitation are presented in Figure 34; with 
forward pitch excitation in Figure 35; and with 
forward lateral excitation in Figure 36. The 
response scale is in ±g per pound of force. 

Vertical vibration prediction from fr'-ward rotor 
vertical excitation in Figure 34 shjws » airly good 
absolute magnitude correlation with test at the 
important 3/rev and 6/rev forcing frequencies 
There is generally an analytical response which 
can be associated with the major test peaks and 
usually the minor ones as well. In the coupled 
direction, i.e. longitudinal motion under vertical 
excitation, the absolute magnitudes, which are 
usually smaller than in the prime directions, are 
reasonably well produced. 


near 17 Hz is reproduced. Vertical motion from 
pitch excitation is acceptable on an absolute basis 
at 3/rev and 6/rev, but the magnitudes of the 
peaks disagree. 

The analytical peak at 32.7 Hz is generally over- 
predicted in amplitude. This implies that the 

proper choice of damping, rather than the constant 
2.5% structural critical damping assumed, would 
improve the correlations. 

Results of the forward rotor lateral excitation are 
in Figure 36. Again, the absolute magnitudes 
are reasonable. On the negative side, the lateral 
peak near 21 Hz is over predicted. Again the use 
of non-constant structural damping would improve 
this situation. 


On the negative side, the very prominent cockpit 
i>ta 52 test response at 28 Hz in the vertical 
direction has no strong analytical counterpart. 

Results of the forward rotor pitch excitation are in 
Figure 35. Comparison of test and analysis here 
gives generally good agreement. Again absolute 
magnitude predictions are good, especially at 3/rev 
and 6/rev. Longitudinal motion at the forward 

hub shows the strong peak near 10 Hz that is 
close to the test peak. Even the secondary peak 


Figure 37 is a bar chart comparing analytical and 
test frequencies. In the cluster of modes from 6 
to 8 Hz, there is one more analytical than test 
mode. Since this analytical mode is an out-of- 
phase engine to engine yaw motion, it may exist 
but be masked within the adjacent aircraft longi- 
tudinal mode at 7.2 Hz. In the cluster of modes 
from 10 to 20 Hz, there is an analytical mode 
corresponding to every test mode. The frequency 
error ranges from near zero to 0.8 Hz for the test 
mode at 11.7 Hz. Above 20 Hz there are more 
analytical than test modes. 
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ANALYTICAL NATURAL FREQUENCIES 
FORWARD VERTICAL/PITCH AND LATERAL TEST CONFIGURATION 
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FiguAt 37. CompaAUon otf Tt&t and Analytical Hatu/ial F/t ijqucncit& 


Correlation Improvements 

A number of items arose from the modeling and 

correlation experience which have the potential for 

further improvement of correlation. 

1. Correct modeling of damping is a major need. 
The current use of a constant assumed value 
of structural damping is not adequate. Some 
form of nonuniformly distributed damping is 
required. 

2. Stringer area is not included in shear area of 
the cross-section, since the usual assumption 
of skin areas carrying all shears is made. 
When summed the shear area of stringers is 
as much as 50% of the skin area. 

3. The upper portion of the splice joints is in 
compression under tg loading and uncon- 
nected stringers may be axially effective. 

4. More thorough modeling of the forward trans- 
mission cover, shaft, bearings and bearing 
clearances may be necessary to obtain a still 
closer match of the mode near 3/rev. 

5. The hub test fixture should be remodeled to 
better reflect elastic effects at the interface 
with the rotor shaft 

6. Masses are distributed to approximately 10% of 
the structural grid points. A finer mesh may 
be necessary to improve higher mode 
predictions. 


A preliminary effort to evaluate some of these 
improvements has been conducted. In Figure 38, 
damping has been adjusted in an attempt to im- 
prove the forced response correlation. Instead of 
using a constant 2.5% structural damping, the 
damping has been varied by mode as indicated in 
the tabulation. 
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Figure. 5f. mtcA oi Modal Vamping on FoAced 
R eAponi* CowitlatLon 


The damping was varied here to obtain the best 
match at the bottom of the response, away from 
the resonance points. 

A second improvement item has been explored. 
Table 2 summarizes the results of a number of ex- 
ploratory runs to Investigate the effect of splice 
joint continuity and stringer shear area. For 
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Table, 1 . E&iect oi S ptice Joint Continuity and 
StAingeA SheaA Anea on Hatu/vat 
TAequency 
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expediency, the stringer shear area was simulated 
by modifying the shear modulus so as to effec- 
tively increase the shear area. The thrust of the 
effort was to raise the baseline analytical fre- 
quency at 10.85 Hz to the test value at 11.7 Hz. 
Tfie chart shows that with all the stringers con- 
tinuous at Stations 160 and 440, the frequency did 
increase from 10.85 to 11.31 Hz. This change in 

splice joint continuity has remarkably little effect 
on the frequency of the remaining modes. 

Next, to represent the actual stringer shear area, 
the shear modulus is increased by a factor of 1.5, 
the frequency of this mode increased to 11.68 Hz, 
almost exactly the 11.7 Hz test value. 

— — MEASURED 

CONTINUOUS SPLICE JOINT 

BASELINE ANALYSIS 

RESPONSE: 5TA, 52 R/H LONG. 

k IQ- jo EXCITATION: FVO. HUB PITCH 



— — MEASURED 

CONTINUOUS SPLICE JOINT AND STRINGER REACTION 


RESPONSE: STA. 52 R/H LONG. 
EXCITATION: FWD. HU8 PITCH 



TiguAt 40. Combined Eiiect oi Splice Joint 
Continuity and StAingeA Shea A 
Reaction 


Forced response runs were then made with these 
two improvements. As representative, look at 
cockpit longitudinal response under forward rotor 
pitching moment excitation shown in Figures 39 and 
40. 

The inclusion of these two, somewhat secondary 
effects, thus have an impressive effect on 
improvement of the correlation. 


Industry Critique of Test and Correlation 

Rather than a series of on site briefings, the 
presentation and critique of the test and corre- 
lation activity was made at a joint meeting of 
industry representatives. The analytical approach 
of modeling the shakers and support systems in 
addition to the basic airframe received favorable 
comments from all attendees. Reasons cited in- 
cluded: (1) verification of normally accepted 

suspension concept, (2) insured one-to-one com- 
parison, and (3) directly addressed interaction 
issue. Overall, the correlation below 20 H2 was 
deemed good. However, the consensus of opinion 
was that the higher frequency range needed more 
work. A finer mass breakdown was considered to 
be a key aspect in improving high freque~' 
correlation. 

Several comments unrelated to specific test results 
are also worthy of mention. One observer sug- 
ges f ed that study of the available results might 
provide guidelines for a realistic validation 

criteria. A second noted that a stronger manage- 
ment commitment to adequate shake testing and 
correlation was needed. 

Ref. 4 reported the details of the ground shake 
test and the correlation effort. Ref. 5 is an 
overall program summary. 
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Conclusions 

• Guides prepared during the planning phase 
enabled proper p!anning / scheduling and 
control of the present modeling effort. 

• Error free demonstration runs for the 
resulting static and vibrations models dis- 
played rational internal loads and reasonable 
natural frequencies and mode shapes. 

• Management enforced cooperation of Design- 
Stress-Weights-Dynamics is key to achieving 
an FEM suitable for internal loads, structural 
member sizing and vibration analysis. 

• Cost of the total effort is 4,430 man-hours or 
5%, 4% is already usual for internal loads; the 
vibration model is another 1%. 

• Satisfactory procedures were developed for 
analysis of the suspended aircraft. Com- 
parison of free and suspended configuration 
indicates only minor differences. 

• Reasonable correlation was obtained between 
test and analytical results. Adequate 
modeling of damping appears as a major 
stumbling block to improved correlation. 

• Nonlinear effects result in test scatter of 
peak responses about the natural frequencies. 
Force level was identified as one source of 
nonlinearity. 

• Significantly improved correlation appears 
possible by including secondary effects such 
as stringer shear area and effective splice 
joint stringer continuity due to Ig loading. 
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Charlie Fredrickson, Sikorsky Aircraft : Dick, I think that was a really neat paper, well 

thought out, and very nicely presented and so on. I’d like to ask you three questions. Did you 
previously do a NASTRAN model on the Chinook, in other words had you previously taken an earlier 
FEM model on the Chinook and upgraded it 'to the D before it actually flew and if you did, how 
did that compare with this well-planned FEM model that you did later on? 

Gabel: Well, we did do that, Charlie, many, many years ago for the A model of the Chinook. 

There was a model built and there was even hardware made to try to tune it and the fact that the 
Chinooks had vibration troubles through their whole life means that it didn't work very well. 

Fredrickson : I was trying to find out if you had upgraded that model for the D. 

Gabel : No, we were too stupid to do that. We actually started from scratch and used the design 

drawings for the D. Since many of the people who did the early model were long gone anyway, it 
didn't really matter. 

Fredrickson: In the actual shake test, how did you actually identify what you considered to be 

natural modes of the aircraft? 

Gabel: Combinations of things: we used the peak, the forced amplitude, we used the 90° phase 

between the amplitudes and the shaker force, we used the frequency circle diagrams--about 3 or 
4 different ways. 

Fredrickson : Okay. I know in my own experience, if you use one or another method and don't use 

kind of a combination, you're liable to miss a few modes along the way. 

Gabel : But then they're not pure because we were shaking with one shaker at a time at one rotor 
head and to get a pure mode you have to have distributed shakers which nobody does any more. 

Fredrickson : Another question about how the shake test was done. Was that a swept sine or 

random input or just exactly what was the methodology behind the shake test itself? 

Gabel : It was a slow sweeping sine. 

Wayne Johnson, NASA Ames Research Center : With the coming switch to composite airframes, do you 

think that's going to make this Job harder or easier? 

Gabel : Different. So far the elements being used are really the same as the stress people have 

been using for the metal elements. They are not going into it layer by layer because of the 
magnitude of the structure. 

Johnson : Do you think the composite structures will have more or less small scale variations? 

It seemed that one of the things you were saying is that small scale variations which are not 
modeled are almost certainly a cause of some of the discrepancies. Do you think composites will 
have more or less of that? 

Gabel: It's hard to say. I would think they might have more because the way they're laid up-- 
it's not quite the same as a rolled out metal sheet. There may be variations in thickness and 
such things that may be more complex. I .night comment that Langley is underway on a continua- 
tion of this program, where Sikorsky, Hughes, and Bell are analyzing their production metal 
aircraft. Since we have already done the first metal one, we're underway on the first composite 
aircraft. We are modeling It and we're going through the same process that's shown hert. 

Bob Wood, Hughes Helicopters : Dick, I'd like to compliment you on a fine presentation. As 

Hughes is one of the participants in it, I Just wanted to bring out one of the values among many 
of the values I think we're finding from this NASTRAN analysis. In the case of the Apache, the 
second vertical bending mode came out to be practically right on our N per rev and we thought 
that was really the problem. But using our NASTRAN [mooel] and taking the percentage of modal 
contributions, it turned out that for the forced response, the primary contributor to the pilot 
and cockpit vibration was a w 'ng-symmetric mode down at 14 Hertz. I think this is one of the 
values we can get out of NASTRAN. 





